Statistical Theory for Incoherent Light Propagation in Nonlinear Media 
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A novel, statistical approach based on the Wigner transform is proposed for the description 
of partially incoherent optical wave dynamics in nonlinear media. An evolution equation for the 
Wigner transform is derived from a nonlinear Schrodinger equation with arbitrary nonlinearity. 
It is shown that random phase fluctuations of an incoherent plane wave lead to a Landau-like 
damping effect, which can stabilize the modulational instability. In the limit of the geometrical 
optics approximation, incoherent, localized, and stationary wave-fields are shown to exist for a wide 
class of nonlinear media. 
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It is well known that the propagation of electromag- 
netic waves and beams in dispersive nonlinear media is 
associated with the phenomena of modulational insta- 
bility and the formation of stationary wave structures, 
i.e., temporal and spatial envelope solitons, cf. |Q. The 
conventional treatment of these phenomena considers co- 
herent wave structures. However, a recent series of exper- 
imental 1 3] and theoretical papers has demon- 
strated the existence of incoherent solitons. In particular, 
spatially incoherent optical solitons have attracted much 
attention. Until now, three approaches have been devel- 
oped to describe partially incoherent wave propagation 
in nonlinear media: (i) the propagation equation for the 
mutual coherence function |3|], (ii) the coherent density 
method |4j] , and (iii) the self-consistent multimode theory 
||. In fact, it has been shown in Ref. || that these three 
approaches are formally equivalent. 

In the present article, we propose a novel, general, sta- 
tistical theory for describing the dynamics of partially 
incoherent optical waves and beams in dispersive and 
nonlinear media. The approach is based on the Wigner 
transform method 0, which was introduced in statisti- 
cal quantum mechanics to describe the dynamics of the 
quantum state of a system in the classical space lan- 
guage. A similar approach has been successfully applied 
in plasma physics in connection with the theory of weak 
plasma turbulence ||, as well as with the description of 
electromagnetic wave propagation in a nonstationary, in- 
homogenous, and relativistic plasma Recently, the 
Wigner function method has been used to analyze the 
longitudinal dynamics of charged-particle beams in accel- 
erators jlO| , and to study the dynamics of Bose-Einstein 



condensates in the presence of a chaotic external poten- 



tial y_L 
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The scope and outline of this article is as follows. 
Starting from the nonlinear Schrodinger equation de- 
scribing the evolution of the slowly varying wave ampli- 
tude in a dispersive medium with an arbitrary nonlinear 
response, we derive the Wigner-Moyal equation for the 
Wigner transform including the Klimontovich statistical 
average. This equation reduces, in the geometrical op- 
tics approximation, to the classical Liouville or Vlasov- 
like equation describing the conservation of optical quasi- 
particles in phase-space, cf. Ref. E3]. To illustrate the 
usefulness of this approach, we consider an application 
to the case of one dimensional (ID) propagation of par- 
tially incoherent light in a nonlinear Kerr medium, and 
investigate the stability of a constant amplitude plane 
wave against small harmonic perturbations. It is found 
that, in addition to the classical modulational instabil- 
ity (MI), a new linear Landau-like damping effect arises. 
This damping is due to the broadening of the Wigner 
spectrum, which is associated with the random phase 
fluctuations. Consequently, it is found that the partially 
incoherent character of the light may suppress the mod- 
ulational instability, in agreement with the result of Ref. 
[PH 14 . Finally, we illustrate an analogue to the BGK 
waves in plasma physics. In the limit of the geometrical 
optics approximation, we derive the Wigner distribution 
functions of a new class of stationary, self-trapped, and 
incoherent wave pulse structures, which may exist in a 
wide class of nonlinear media. 

The approach proposed here sheds new light on the 
physics behind the recently reported results in this 
field. It also allows the formulation of conceptually new 
problems regarding the dynamics of partially incoherent 
waves and beams in nonlinear media. In addition, us- 
ing the Wigner transform for studying this kind of phe- 
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nomena establishes new connections to quantum mechan- 
ics, plasma physics, mesoscopic physics, image and signal 
processing, and mathematics. Furthermore, the Wigner 
transform, apart from serving as a very convenient math- 
ematical tool, is widely used in its own right in the field 
of optics, as a supplement to the envelope wave function 

As our starting point, we assume that the three- 
dimensional (3D) wave propagation in a dispersive (or 
diffractive) nonlinear medium is described by a system 
of coupled model equations for the slowly varying com- 
plex wave amplitude ip(t, r) and the nonlinear response 
function of the medium n(t,r): 

i i§t +Vg ' V ^) ^ + 2 V ^ + ntf; = 

T m —+n = KG((ip*ij)) 
at 

For convenience, we use t and r as the evolution and 
spatial dispersive variables, respectively; v g is the group 
velocity, f3 is the diffraction or second order dispersion 
coefficient, k is the nonlinear coefficient, and the func- 
tion G({if)*ijj)) characterizes the nonlinear properties of 
the medium. The bracket (. . .) denotes the statistical en- 
samble average. The relaxation time of the medium re- 
sponse function, r m , is assumed to be much longer than 
the characteristic time of the statistical wave intensity 
fluctuations, r s . Assuming also that r p 3> r m , where t p 



is the characteristic time scale of the (deterministic) wave 
amplitude variation, we can approximate the medium re- 
sponse function as n ~ kG ({ip*ip)). The system (Q) then 
reduces to a generalized nonlinear Schrodinger equation: 



i^ + |V 2 ^ + «G((^))^ = 0, (2) 

where the coordinate system has been transformed to the 
reference system moving with the group velocity v 9 , and 
r is the reduced distance according to r — > r — v g t. 

We stress that the independent variables in Eq. (^) 
should be adopted to suit the particular situation. For 
instance, to describe the 2D CW spatial optical solitons 
one would use z as the evolution variable and restrict the 
nabla operator to the transverse dimensions, while an 
equation similar to Eq. (|J) with z as evolution variable 
and d 2 /dt 2 instead of the nabla operator would describe 
temporal ID solitons. 

Replacing id/dt by oj and — iV by p yields the corre- 
sponding nonlinear dispersion relation: 



co=^p 2 -KG((r^)- (3) 

We will now apply the Wigner transform method 
0. The s-dimensional Wigner transform, including the 
Klimontovich statistical average, is defined as 



1 f +OQ 

p(p, t, r) = — — / d s £ e*P-« (r (r + £/2, t) ifi (r - £/2, t)> , 



(4) 
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which satisfies (r/>* (r, t) (r, t)) = f_°° d s pp (p, t, r). 
Applying this transform, with dimensionality s = 3, to 
Eq. (g), we obtain the following Wigner-Moyal equation 
for the evolution of the Wigner distribution function: 




0. (5) 



The arrows in the sine differential operator indicate that 
the derivatives act to the left and right, respectively. The 
sine operator is defined by its Taylor expansion: 
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For those situations where the geometrical optics ap- 
proximation is valid, we can neglect the higher order 
derivatives in Eq. (ra). This approximation is valid for 



Ap • Ar 3> 2-7T, where Ap is the local width of the Wigner 
spectrum and Ar is the width of the medium response 
function n(t, r). We emphasize that this is a rather strong 
limitation, which can only be valid in the long-wavelength 
limit. Retaining only the first term in the expansion (o), 
we obtain a Vlasov-like equation: 



dp a dp dG((w 
at or or 



op 



(7) 



Equation (Q) implies the conservation of the number of 
optical quasi-particles in phase-space. Using the Liou- 
ville theorem we recover the canonical Hamilton equa- 
tions of motion for the quasi-particle with mass /3, or, 
equivalently, the ray equations of the geometrical optics 
approximation, where r and p are the canonical variables 
and u> is the Hamilton equation defined by Eq. (|^): 
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¥ = 

du> 

dr 



0P, 



dG{(r^)) 
dr 



(8) 



Note that Eg. ( H) is similar to the radiation transfer equa- 
tion used in [|^|l7). 

As an example of the applicability of the Wigner 
transform approach, we consider the MI of a ID plane 
wave with a constant amplitude in a nonlinear Kerr-likc 

medium, for which G (J^\ip\ 2 X\ = {^IV'I 2 ^- F° r the case 
of coherent light propagating in a ID Kerr medium, it 
is well known from conventional stability analysis that 
a perturbation of the monochromatic stationary solution 
ip(x,t) = ipo exp (inipQij experiences a modulational in- 
stability when (3k > and K 2 < AkiPq/(3, where K is the 
wave number of the perturbation. The instability growth 
rate is given by 
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(9) 



To investigate the effect of the incoherence, we assume 
a Wigner distribution function of the form p(jp, t, x) = 
Po(p) + Pi ex P [i {Kx — fit)], with p > \pi\- Here, 
Pa{p) is the background distribution function correspond- 
ing to the plane wave with a complex amplitude ^ = 
1(>0 exp [inipfit + i<f>(x)] , implying f_™p (p)dp = Vo- 
The incoherence is modeled by the randomly varying 
phase term <fi(x). Studying the linear evolution of the 
perturbation p\, we obtain from the linearized Wigner- 
Moyal equation (||) the following dispersion relation: 



Wigner spectrum because of the interaction between dif- 
ferent parts of the spectrum. This spectral redistribution 
counteracts the MI. Similar phenomena occur in connec- 
tion with nonlinear propagation of electron plasma waves 
interacting with intense electromagnetic radiation j^, ^| , 
nonlinear interaction between random phase photons and 
sound waves in electron-positron plasmas Jl7t , and the 
longitudinal dynamics of charged-particle beams in ac- 
celerators [ [To[ . 

It is interesting to note that for a coherent wave, i.e., 
a delta-shaped background distribution pa(p) = ^qS(p), 
and for /3k > 0, Eq. (0) gives exactly the modula- 
tional instability growth rate defined b y E q. (Q). For 
the same case, the dispersion relation ( |ll| ) yields fl — 

i (Pk) 1 ^ 2 tpoK , which is identical to the growth rate of 
the MI as given by the expression (0) in the limit of long 
wavelengths, i.e., when K 2 <C 4k?/>5//3. 

To illustrate the incoherent case, we assume that <p{x) 
is described by the following autocorrelation function: 



(exp [-i(f> (x + y/2) +i<j)(x- y/2)]) = exp (-p \y\) , 

(12) 

where p is the correlation length. The corresponding 
Wigner function has a Lorentzian shape: 



Po(p) 



7T p Z +Pq 



(13) 



and the dispersion relation (|To|) yields the following exact 
result for iKip 2 /{3K 2 > 1: 



fl _ K (Aki^I 
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(14) 



po (p + K/2) -p (p~ K/2) 
K(p-fl//3K) 



dp = 0. (10) 



The corresponding dispersion relation following from the 
linearized Vlasov-like equation (Q) has the form 



1 



dpg/dp 
(p - fl/pK) 



dp = Q 7 



(11) 



which can also be directly obtained from Eq. ( |l0| ) by tak- 
ing the limit of small K. The relation ([ll]) is similar to 
the dispersion relation for electron plasma waves, which 
is well known to contain the effect of the Landau damp- 
ing. In general, the kinetic integrals in ( |To| ) and jll| ) 
can be represented as the sum of a principal value and a 
residue contribution, where the latter leads to a Landau- 
like damping of the perturbation. This stabilizing effect 
is not an ordinary, dissipative damping. It is rather an 
energy-conserving self-action effect within a partially in- 
coherent wave field, which causes a redistribution of the 



which is similar to the result obtained in Rcf. p3| . 

Equation ( |l4| ) clearly shows the stabilizing effect of the 
Landau-like damping due to the finite width pq of the 
Lorentzian spectrum (or the finite correlation length of 
the wave phase). In fact, if the width of the Lorentzian 
spectrum, po, satisfies the relation po > p c = K c /2, 

1/2 

where K c = (4k^o//3) is the cut-off wavelength of 
the MI, the instability is completely suppressed for all 
wave numbers K of the perturbation. In other words, if 
the Landau-like damping, induced by the broadening of 
the Wigner spectrum due to the partial incoherence of 
the wave, is strong enough, it can overcome the coherent 
growth associated with the MI. This result has also been 
verified experimentally; see, e.g., Ref. Jl4| . 

Finally, we will demonstrate the existence of a class of 
self-trapped and partially incoherent wave pulse solutions 
to the stationary Vlasov-like equation 
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dp 
dx 
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dp 



(15) 
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for an arbitrary nonlinear function G (^l^| 2 ^) • Accord- 
ing, to the Jeans theorem, c.f. |^| , the solution of Eq. 
( |l5| ) can be expressed as an arbitrary function of the 

Hamiltonian H = j3p 2 /2 — kG ((lV-'| 2 ^) • Thus, we have 
p = p s (H), where p s is an arbitrary function. The quasi- 
particles are trapped in the nonlinear potential for —7 < 

1 /2 

p < 7, where 7 = (2kG ((|V>| 2 )) / P\ and k/3 > 0, 

and consequently the condition dp po(p) = (\tp\ 2 
leads to an integral equation having the following solu- 
tion for p s (c.f. |lql ): 



1/2 



dH' dF (H') 
(-H-H') 1/2 dH> 



> (16) 



where F (G) = G^ 1 (9) with 6 = G ((\ip\ 

(/3p 2 /2 — H^j I k. For instance, for the Kerr nonlinear- 
ity we have F (8) = 9 and Eq. @ yields p s (H) = 
(2/irn) y/—/3H/2. Equation ( [l6| ) describes the Wigner 
distributions of quasi-particles trapped in a collectively 
produced, stationary, partially incoherent, and localized 
wave structure. These structures are similar to the large- 
amplitude BGK waves in a plasma; see p8[ . 

In conclusion, we have proposed a general, statisti- 
cal approach for the theoretical analysis of the propa- 



gation of partially incoherent optical waves and beams 
in dispersive media with arbitrary nonlinearities. The 
approach is based on the Wigner transform method in- 
cluding the Klimontovich statistical average. The derived 
Wigner-Moyal equation determining the evolution of the 
Wigner distribution function represents a generalization 
of the Vlasov-like equation, which is only valid within the 
geometrical optics approximation. The Wigner-Moyal 
equation clearly shows that the number of optical quasi- 
particles is not conserved in phase-space beyond the va- 
lidity of the geometrical optics approximation. Using 
the Wigner-Moyal equation, we have carried out a linear 
stability analysis for small perturbations on a constant, 
ID, partially incoherent background in a nonlinear Kerr 
medium. The theory reproduces the exact expression for 
the MI growth rate of a coherent wave, but also includes 
a linear Landau-like damping effect associated with the 
broadening of the Wigner spectrum due to partial wave 
incoherence. This damping effect explains the previously 
reported incoherent suppression of the modulational in- 
stability. 



Acknowledgments 

The authors are grateful to D. N. Christodoulides, M. 
Segev and A. Hasegawa for valuable discussions. 



[1] V. I. Karpman and E. M. Krushkal, Sov. Phys. JETP 
28, 277 (1969); T. Tanuti and H. Washimi, Phys. Rev. 
Lett. 22, 454 (1969) 

[2] M. Mitchell et al, Phys. Rev. Lett. 77, 490 (1996); 
M. Mitchell and M. Segev, Nature (London) 387, 880 
(1997); Z. Chen et al, Science 280, 889 (1998) 

[3] G. A. Pasmanik, Zh. Eksp. Teor. Phys. 66, 490 (1974) 
[Sov. Phys. JETP 39, 234 (1974)]; V. V. Shkunov and 
D. Z. Anderson, Phys. Rev. Lett. 81, 2683 (1998) 

[4] D. N. Christodoulides et al, Phys. Rev. Lett. 78, 646 
(1997); D. N. Christodoulides et al, Opt. Lett. 22, 1080 

(1997) ; T. H. Coskun et al, Opt. Lett. 23, 418 (1998) 
[5] A. Hasegawa, Phys. Fluids 18, 77 (1975); 

M. Mitchell et al, Phys. Rev. Lett. 79, 4990 (1997); 
D. N. Christodoulides et al, Phys. Rev. Lett. 80, 2310 

(1998) ; M. I. Carvalho et al, Phys. Rev. E 59, 1193 

(1999) ; T. H. Coskun et al, Phys. Rev. Lett. 84, 2374 

(2000) ; N. Akhmediev, W. Krolikowski, and A. W. Sny- 
der, Phys. Rev. Lett. 81, 4632 (1998); W. Krolikowski, 
N. Akhmediev, and B. Luther-Davies, Phys. Rev. E 59, 
4654 (1999); A. A. Sukhorukov and N. N. Akhmediev, 
Phys. Rev. Lett. 83, 4736 (1999) 

[6] D. N. Christodoulides, E. D. Eugenieva, T. H. Coskun, 
M. Segev, and M. Mitchell Phys. Rev. E 63, 035601 

(2001) . 

[7] E. Wigner, Phys. Rev. 40, 749 (1932); J. E. Moyal, Proc. 



Cambridge Philos. Soc. 45, 99 (1949); see also the text- 
books in statistical quantum mechanics 
[8] A. A. Venedov, Theory of Turbulent Plasma, Moskva 
1965, Israel Program for Scientific Translations, 
Jerusalem 1966 
[9] N. L. Tsintsadze and J. T. Mendonca, Phys. Plasmas 5, 
3609 (1998); Phys. Rev. E 62, 4276 (2000) 
[10] D. Anderson et al, Phys. Lett. A 258, 244 (1999); 
R. Fedele, D. Anderson, and M. Lisak, Physica Scripta 
T84, 27 (2000) 
[11] S. A. Gardiner et al, Phys. Rev. A 62, 023612 (2000) 
[12] A. I. Akhiezer, I. A. Akhiezer, R. V. Polovin, 
A. G. Sitenko, and K. N. Stepanov, Plasma Electrody- 
namics, Pergamon Press, 1975 
[13] M. Soljacic et. al, Phys. Rev. Lett. 84, 467 (2000) 
[14] C. Anastassiou et. al, Phys. Rev. Lett. 85, 4888 (2000) 
[15] For a long list of references, see the JOSAA feature issue: 
Wigner Distributions and Phase Space in Optics, edited 
by G. W. Forbes, V. I. Man'ko, H. M. Ozaktas, R. Simon, 
and K. B. Wolf, JOSAA 17, 2274 (2000) 
[16] R. Bingham, J. T. Mendonca, and J. M. Dawson, Phys 

Rev. Lett. 78, 247 (1997) 
[17] P. K. Shukla and L. Stenflo, Phys. Plasmas 5, 1554 (1998) 
[18] I. B. Bernstein, J. M. Green, and M. D. Kruskal, Phys 
Rev. 108, 546 (1957). 



